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Abstract 

Using the Sakai-Sugimoto model we study the effect of an external magnetic field 
on the dynamics of fundamental flavours in both the confined and deconfined phases 
of a large N c gauge theory. We find that an external magnetic field promotes chiral 
symmetry breaking, consistent with the "magnetic catalysis" observed in the field the- 
ory literature, and seen in other studies using holographic duals. The external field 
increases the separation between the deconfinement temperature and the chiral sym- 
metry restoring temperature. In the deconfined phase we investigate the temperature- 
magnetic field phase diagram and observe, for example, there exists a maximum critical 
temperature (at which symmetry is restored) for very large magnetic field. We find 
that this and certain other phenomena persist for the Sakai-Sugimoto type models 
with probe branes of diverse dimensions. We comment briefly on the dynamics in the 
presence of an external electric field. 
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1 Introduction 



Since the advent of the AdS/CFT correspondence [H [5], [3], (see a review in ref.[l]) there 
has been considerable refinement of the methods for studying strongly coupled large N c 
gauge theories. Much effort has been spent studying such systems at both zero and finite 
temperature, constructing specific "holographic" models intended to capture key features 
of QCD at strong coupling, such as the confinement /deconfinement phase transition, chiral 
symmetry breaking, and possible novel phases that may be of relevance to experiment and 
observation. 

The Sakai-Sugimoto model, as described in ref. [7] is one such construction which cleanly 
realizes chiral symmetry breaking and deconfinement. The supergravity background of 
this model is constructed of near- horizon geometry of N c ZM-branes, following ref. [5]. The 
study of Nf flavour D8 branes in this background when Nf <C N c reveals a nice geomet- 
ric realization of chiral symmetry breaking. The flavour branes do not backreact on the 
background geometry in this probe limit and therefore studying their dynamics using the 
Dirac-Born-Infeld (DBI) action (including a Wess-Zumino term, if necessary) suffices to 
capture the corresponding gauge theory dynamics of fundamental flavours in an analogue 
of the quenched approximation. 

In this (DBI) regime it is possible to capture general gauge theory features such as the 
phase diagram for temperature vs chemical potential by considering probe brane in finite 
temperature supergravity background and exciting specific gauge field on the world-volume 
of the probe brane. Previous such studies including the non-zero chemical potential in this 
model have been carried out in e.g., refs. [T2l [T^ [T5 | [T6] . Here, we will introduce an external 
magnetic and electric field. 

A clear method for introducing a background magnetic field has been previously discussed 
in the D3/D7 model in ref. [26J. The authors consider pure gauge 5-field in the supergravity 
background, which is equivalent to exciting a gauge field on the world-volume of the flavour 
branes corresponding to a magnetic field. 

We find that the presence of magnetic field promotes the spontaneous breaking of chiral 
symmetry. This is expected from the field theory perspective and is widely recognized as a 
sort of "magnetic catalysis" for chiral symmetry breaking (see e.g., ref. [30]). Further study 
of the phase structure of this model reveals the existence of a finite critical temperature 
(for restoration of chiral symmetry) for large magnetic field. We analyze a number of other 
physical quantities such as the latent heat and relative magnetisation associated to the phase 
transition. 
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We also find that our phase structure is rather generic for the Sakai-Sugimoto type holo- 
graphic models where we consider the dynamics of probe Dp-brane in Z)4-brane background. 
We address some of the physics of an external electric field in the model. We find that in 
the symmetry-restored phase an external electric field drives a current in the gauge theory 
due to pair creation, and the symmetry-broken phase does not conduct. However we have 
not considered the presence of baryons in our set-up, which could give rise to a non-zero 
current in the phase where chiral symmetry is broken. 

This paper, is organised as follows: We briefly review the Sakai-Sugimoto construction 
in section 2. In section 3 we perform the analysis of probe D8-branes in the presence of 
magnetic field, while in section 4 we discuss similar results for general probe Dp-brane. In 
section 5 we comment on some aspects of the physics of an external electric field, concluding 
in section 6. 

Note added: When this paper was being prepared we became aware of ref. [28], in which 
authors have studied related physics. 



2 The Sakai-Sugimoto Construction 



The Sakai-Sugimoto model[7] consists of near-horizon limit of N c Z)4-branes wrapped on 
a circle of radius R in the x 4 direction with anti-periodic boundary condition for fermions. 
The D 4- branes are intersected in the compact x 4 direction by Nf D8-branes at x A = — ~ 
and Nf D8-branes at x 4 = |f (with the constraint that L < irR). This is dual to a 
(4 + l)-dimensional SU(N C ) Yang-Mills theory with gauge coupling constant g 5 ; left and 
right handed quarks are introduced by the flavour D8 and D8-branes in the probe limit 
that share three spatial directions with the ZM-branes. The flavour branes introduce a 
global flavour symmetry U(Nf)i x U(Nf)n as seen from the (4 + l)-dimensional ZM-brane 
worldvolume gauge theory. This global symmetry is identified with the chiral symmetry 
(non-abelian) of the effective (3 + l)-dimensional gauge theory with chiral fermions. In 
the probe limit, the dynamics of the flavour branes is described by the DBI action in the 
background of N c D4-brane geometry. The background metric of D4-brane is obtained from 
type IIA supergravity and is given by 

ds2 = (£ I ) 3/2 (-* 2 + ^' + /W(^)V (£f ^(f) • 



RnJ ' " Va ^ ' \ U 
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Here x l are the flat 3-directions, t is the time coordinate, x 4 is the spatial compact circle, 
Q4 are the S 4 directions and u is the radial direction. l s is the string length, g s is the string 
coupling; V4 and €4 are the volume and volume form of S 4 respectively. Also, <fi is the dilaton 
and Fu) is the RR four-form field strength. To avoid a conical singularity in the {x 4 ,u} 
plane one should make periodic identification: 



4tt / R 



?3 \ 1/2 



This endows the background with a smooth cigar geometry in the {x 4 , u} plane. The 
radial parameter R? m is given by 

R 3 m = ng s N c ll = irXa' , (3) 

where A is the 't-Hooft coupling. This construction has three dimensionful parameters 
#5, L and R, where the L and R have been defined above. The five- dimensional gauge 
coupling is given by g\ = (2ii) 2 g s l s . The four dimensional gauge coupling can be obtained 
by dimensional reduction yielding g\ = glJI-nR. The five dimensional 'tHooft coupling is 
defined to be A = (gfiV c )/47r. The gravity picture is valid for small curvature which amounts 
to A R, namely for strong 'tHooft coupling. It should be noted that due to the presence of 
a varying dilaton in equation (CQ) the type IIA supergravity background becomes unreliable 
in the far UV and we need to lift it to M-theory for a possible UV-completion. 

To consider the finite temperature version of the model we need to Euclideanise the 
background in equation ([!]). This can be achieved by compactifying the time direction, t, 
on a circle and identifying the period with inverse temperature (3. In this case, the x 4 circle 
shrinks away at u = Ut but the t circle is fixed. One can also construct a finite temperature 
version by interchanging the role of the t and x 4 circles so that now time circle shrinks away 
at some value u = Ut but the x 4 circle remains fixed. It is easy to see that both these 
constructions have the same asymptotic behaviour. These are the only known Euclidean 
continuations of the background in equation (0Q) with the right asymptote. It is known (e.g, 
in ref. [9]) that for T < 1/2ttR (i.e low temperature), the background with the x 4 circle 
shrinking dominates, where for T > 1/2ttR (i.e., high temperature) the background with 
the t circle shrinking dominates; and this geometric transition between the two background 
corresponds to the confinement /deconfinement transition. 

So for low temperature the relevant background is given by equation ([1]) with the time 
coordinate periodically identified with period (3 = 1/T, where T is the temperature. The 
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high temperature background is given by 

ds 2 = (!) 3/2 {dxtdx* + f(u)dt 2 + (dx 4 ) 2 ) + (|)' 3/2 + u 2 dti^ , 

All the parameters are given by the same formula as equation ([1]). The dilaton, RR 4-form, 
_Rd4 are also given by the same formula as equation ([1]). Now one can introduce the flavour 
brane-anti-brane system in the probe limit, namely Nf <C N c . In this limit the probe branes 
do not backreact on the geometry and the classical profile of the probe is solely determined 
by the Dirac-Born-Infeld action. We will consider the following ansatz for the flavour D8-D8 
branes: 

{t, Xi, x A = r, O4, u = u(t)} . (5) 

For notational convenience we rename x A coordinate to be r, and we note that the co- 
ordinates in the parenthesis should be understood as the worldvolume coordinates of the 
D8/D8-brane. 

3 The Probe Brane Analysis 

Many aspects of the finite temperature physics have been studied before, e.g. in refs. [9].|llj 
and references therein. To introduce external magnetic field we follow the procedure adopted 
in ref. [2S]- We consider the presence of a pure gauge 5-field given by, B 2 = Hdx 2 A dx 3 . 
As far as the DBI action is concerned it can be easily seen that such a choice is equivalent 
to exciting a gauge field A 3 = Hx 2 on the worldvolume of the probe brane. 

3.1 The Low Temperature Background 

In this case the relevant background is given by equation ([1]). With the probe brane ansatz 
in equation (jHJ) the induced metric on the worldvolume of D8/Z^8-brane is given by 
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where u' = du/dr; also we drop the negative sign in front of the time coordinate as we are 
considering the Euclidean metric. The DBI action is given bjl] 



S 



D8 



/'8 



J d 9 ie-^det(P[G^ + B^}) = C J drC(u,u') 



(7) 



where C = ^VsiV^ / g s T . Recall that we put the flavour branes (D8/D8) with the asymp- 
totic condition that as u — > oo, r — > ±L/2; where the ± corresponds to the D8 and Z)8-brane 
respectively. Also note that the coordinate r is restricted from — tcR to +nR. 

One can immediately note from equation (|7j) that C(u,u f ) is independent of r, therefore 
the Hamiltonian corresponding to r will be a constant of motion. Carrying out the following 
Legendre transformation we should have 

fdC(u, v!) 



It 



du' 



C(u, u') = const. 



So the first integral of motion that follows from equation (jSJ) is given by 



u 



1 + H 2(R^fy f{u 



r/ 4 



1 + H 2 



R 



Dl 



Uo 



(9) 



We have rewritten the constant in the right hand side in a convenient way. Note that 
Uo is the minimum value of u that the probe brane can reach satisfying u'\ u= u = C@. For 
zero background magnetic field this set up reduces to the low temperature case analyzed in 
ref. j9]. Let us now focus on the solution for the probe brane profile. 

We will compare the behaviour of the brane profile in the presence of magnetic field to 
the case when it is turned off. For notational convenience we define the following: 

Ukk 



u 

y = w 



Vkk 



Un 



R 



D4 



U n d 



Ual 



With the above redefinitions we can obtain the difference in slope of the profile in presence 
and in absence of magnetic field as 

H 2 d 3 



.'2 



U H U H=0 



f(y? fiy) 



-3 y 



1 + H 2 d 3 



1 



: It is a straightforward exercise to check that in this case the Wess-Zumino term does not contribute to 
the probe brane action. 

2 We assume that the brane-anti-brane pair join smoothly, which implies that there is no resultant force 
present at the point where they meet. Typically this would mean that there is no other source (e.g., a 
baryon vertex or a bunch of F-strings) present at this point. 
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with y G [1, oo]. So we get that \u' H \ < \u' H=Q \ for each value of y. This in turn means that 
the magnetic field bends the profile of the D8/D8 brane and therefore forces the brane- 
anti-brane pair to join closer to the boundary (and hence break chiral symmetry) for fixed 
asymptotic separation. 

We can study this explicitly as follows. The brane-anti-brane separation at the boundary 
(u — > oo) is given by 



dr 



2 J Ju 0(H) u 'h \/ u o{h) 



i+jf 2 d a /(i) y 1 



f(y) 

l( H )(y)dy . (10) 

Clearly putting H = we get the corresponding separation when the background magnetic 
field is switched off. 



l R 3 £ r y- 3/2 dy 



2 VUoJi /(y) (^-l) 1/2 ^ 



dyl H =o ■ (11) 



For the same asymptotic separation magnetic field changes the brane profile's point of 
closest approach Uq(h)- We can compare Uq(h) and Uq. Equating equation ( Tl0|) and equa- 
tion (flTT) one gets 

1 Uq{h) = f™l H =ody 
U f™l (H) dy 

Some algebra shows that Th=o > 1{H) for all y, so the ratio on the right hand side is greater 
than or equal to one, which also means that Uq(h) — Uq. Therefore for the same asymptotic 
separation the magnetic field can only help to join the brane-anti-brane pair favouring chiral 
symmetry breaking. This is pictorially represented in figure [TJ 

We can extract more qualitative features in appropriate limits. To do so, let us rewrite 
equation ( TlOl with the change to variable z = y~ 3 . With this equation (TlOl becomes 

L_R 3 £ fZ .,3 7Z , mm f 1 (l-y 3 KK zy l z-^d 



y 3 KK)^ + H 2 d 3 ) / v yKK 1 (13) 

Now small asymptotic separation corresponds to large values of Uq which means yxK *C 1. 
So for small L and weak magnetic field (1/d 3 / 2 ^> H), the leading behaviour of the separation 
is given by (using equation (TlBl ). L ~ R^/^/Uq. This is same as the leading behaviour in 
zero magnetic field case in ref. [9]. However, for strong mag netic field {l/d?' 2 < H), the 
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u 




Figure 1: The dashed U-shaped curve represents a profile in vanishing background field and the 
solid U-shaped curve represents a profile when a non-zero magnetic field is present. These profiles 
are obtained by numerically solving the equation of motion for the probe brane. 



leading behaviour obtained from equation (fT3j) is given by, L ~ R^ da H/Uq. So for fixed value 
of Uq the asymptotic separation scales with the applied magnetic field strength H. This is 
however true only in the Ukk *C 1 limit. The general dependence is more complicated and 



a numerical study yields figures 2(a) and 2(b) 
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(b) 



Figure 2: The dependence of the asymptotic separation between the flavour branes in units 



of R 3 ^/VUo with ukk and magnetic field in the low temperature phase. Figure 2(a) shows 



the dependence of L with ukk for different magnetic field strength; blue (bottom most) curve 
corresponds to H = 0, red (middle) corresponds to H = 1.0, green (top most) corresponds to 
H = 2.0. Figure 2(b) shows the behaviour of L with applied magnetic field for different values of 
Ukk'i blue (bottom most) curve corresponds to ukk = 0, red (middle) corresponds to ukk = 0.5, 
green (top most) corresponds to yxK = 0.7. We have set d = 1. 



From figure 2(a) we can see that in the low temperature phase the asymptotic separation 
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increases as the point of joining (i.e, Uq) of the flavour branes decreases. This means that 
the end points of the D8/D8 move further and further away as we go deeper and deeper in 
the core. The role of magnetic field is to further increase this asymptotic distance for a given 
Ukk- However as we approach y KK = 1, the magnetic field does not affect the separation of 
the flavour branes any more, since all curves start converging rapidly near ukk = 1- This 
is the point where the background geometry ends precisely where the flavour branes join, 
therefore all probes for any asymptotic separation should end at this point irrespective of 
where they start from at infinity. This is consistent with earlier studies in ref. [13J. 



On the other hand, figure 2(b) shows that for a fixed value of Ukk magnetic field can 
increase the asymptotic separation, but not indefinitely. This means that for a fixed Dkk 
as we increase the magnetic field the flavour branes move further away from each other, but 
for high enough magnetic field this separation saturates and becomes insensitive to further 
increment of the magnetic field. We will see later that such saturation shows up in other 
physical quantities also. The role of yKK here is to shift each curve upwards as we increase 
its value. 

In ref. [8] it was noted that the special case of i/kk = (meaning when the radius of 
the spatial circle goes to infinity and the compact direction becomes a flat direction) the 
supergravity background is dual to a non-local NJL model in which the separation scale 
between the brane-anti-brane pair (denoted as L here) determines an effective coupling for 
a four fermi interaction term. As we have seen the magnetic field affects the asymptotic 
separation and therefore tunes the effective coupling. 

Note that in this model the bare quark mass is always zero as there is no separation 
between the flavour and the colour branes at the boundary; however since the branes join at 
some length scale Uq > Ukk in the core one can consider a string stretching from u = Ukk to 
u = Uq. The mass associated with the string can be identified to be the effective constituent 
quark mass as argued in ref. pj. If we denote this mass by M q , then M q = Ju° K y/9u9uu 
, where Uq has to be determined from equation fflBT) for given L. This turns out to be a 
self-consistency equation for Uq. This is analogous to the Gap equation in the field theory 
context (e.g., in ref. p]). The self-consistency equation turns out to be 

Uo = \^ (1 - vIk) {Ul + H 2 R 3 m ) I(U , H) 2 , where 



I(U ,H) 



z 



(z-W (1 - y KK z) (U$ + H*B? m ) - (1 - y KK ) (^o 3 + H 2 R% A )f 2 ' 

(14) 



Now equation fTl4l) can be solved perturbatively, i.e., starting with an initial value for the 
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parameter Uq we can determine the next order approximation to Uq using equation ffl~4l ; 
and we continue until the desired accuracy has been achieved. It is straightforward to guess 
the initial value of Uq. Plugging in H = in equation (jHJ) we should get the constituent 
mass for the low temperature case. This can serve as the initial guess for small magnetic 
fields. Once Uq is known for small magnetic fields, it can be used as the initial guess for 
successively higher values of magnetic fields. Thus we obtain the dependence which is shown 
in figure below. 




H 



0.5 1 1.5 2 



Figure 3: The dependence of constituent quark mass (measured in units of (2ira') 1 ) on the 
external magnetic field in the low temperature phase. We have fixed Ukk = 0.4 and Rd4 = 1. 



It is expected that the mass of the vector and pseudoscalar meson would monotonically 
increase as the constituent quark mass increases. Figure [3] therefore should capture the 
behaviour of massive vector and pseudoscalar meson spectra in presence of magnetic field. 

Note that for this background disjoint brane pair do not exist: The constant r-solutions 
namely, r = —ixR/2 and r = +nR/2 join at u — Ukk- This is because of the cigar geometry 
of the background in the {r, u] submanifold and the fact that r-circle is wrapped by the 
probe branes. The brane pair must join together since there is no place in the geometry for 
them to end separately. So the only configuration possible in the low temperature phase 
breaks chiral symmetry by reducing the global U(Nf) L x U(Nf) R to the diagonal U(Nf). 
This is geometrically understood as the joining of the flavour 8-branes. 

We can compute the action for this brane by substituting for the profile function into 
equation (J7J). The result is 

' i,,N 3 m m) Jo [(1 - y KK z) (1 + zH 2 d s ) - zW (1 - y KK ) (1 + HW)} 1 ' 2 ' 

(15) 
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Note that for small L, this energy is proportional to the vacuum energy of the system. 
Recalling the behaviour of L for small yxK and weak magnetic field we get Sds ~ 2(772^/ L 7 . 
This is similar to the leading order behaviour in absence of magnetic field obtained in ref . [H] . 
On the other hand, for small ukk and strong magnetic field we get Sd 8 ~ 2CR 3 D4 H 1 ^ 2 / L 7 ^ 2 . 
So sufficiently high magnetic field changes how the vacuum energy blows up as L — > 0. The 
general dependence is more involved which we do not pursue here. 



3.2 The High Temperature Background 



Recall that the high temperature background is given by equation (jlj); and we again use the 
same ansatz for the probe given by equation ([5]). As before the profile of the probe brane 
is completely determined by the DBI action, from which a first integral of motion can be 
readily obtained to be 



u 



(l+H 2 (^fffiu) 



H 1 



u 




(16) 



For convenience we use the dimensionless variables defined in equation (Tl0j) along with the 
new variable yx = Ut/Uq. Starting from the first integral of motion in equation (1161) it is 
easy to verify that the finite temperature analogue to equation (fTUj) takes the following form 



.'2 



.'2 



U H U H=0 



y\ 3 f(yf 



H z d 



2j3 



dJ f(l) * 1 + H 2 d 3 \y 



(17) 



This also suggests that \u' H \ < \u' H=0 \, leading us to the same conclusion that the magnetic 
field helps bending the branes. The analogue to equation (TTOT) and (TTT|) now take the 
following forms 



L 

2 

L 
2 



R 3/2 



U 3/2 dy 



0(H) Ji 
R 3/2 



VfW) 



v 7 ^ 



0(H) Jl 



/(i) i+h 2 <p y 1 

y~ z l 2 dy 



VfW) 



/(y)„.s 



/(i) 



1/2 



1/2 



I dy 



l H dy 



leading us to a similar conclusion as the low temperature case. This is pictorially represented 
in figure HI 

As before one can extract the dependence of the asymptotic separation in the small yr 
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T 



Figure 4: The dashed U-shaped curve represents a profile in vanishing background field and the 
solid U-shaped curve represents a profile when a non-zero magnetic field is present. The straight 
(red) solution does not have any qualitative change in presence or absence of the external field. 
These profiles are obtained by numerically solving the equation of motion for the probe brane. 



limit. Again with the change to variable z = y 3 we get 



L 
2 



R 



3/2 
DY 



0(H) 



y T ) (1 + H 2 d 3 ) 



(1 



Vt z ) 



-1/2 



(19) 



- (1 - y*) (1 + H 2 d 3 ) 



So from equation (TT9T) one can see that for small and weak magnetic field L ~ R 3 ^/ ^/Uq(~h)] 



The general de- 



whereas for small yr and strong magnetic field we get L 
pendence has been studied numerically and the result is shown in figure 5(a) and 5(b) 



In figure 5(a) we see that the asymptotic separation decreases as yx increases, consistent 



with studies in ref. p3]. The role of magnetic field is to shift the curves upwards, namely 
to increase the asymptotic separation. However this effect vanishes as yx approaches its 
maximum value and the separation becomes insensitive to the background magnetic field. 



Figure 5(b) shows a similar behaviour as the low temperature case. The separation at the 
boundary becomes higher and higher for increasing magnetic field, but for high enough field 
the flavour branes at the boundary tend to not sense any further increment (therefore a 
saturation is obtained). In this case however, curves for different yx may intersect each 



other (as shown in figure 5(b) ) unlike the low temperature case. 



The joining of the flavour branes inside the core can be associated with the effective 
constituent quark mass. This corresponds to a self-consistency equation for Uq as in the low 
temperature case. The equation can again be solved using the same perturbative approach 
and the results are summarised in the figure El This behaviour of constituent mass is valid 
when the curved solutions are the lowest energy solutions (the chiral symmetry broken 
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Figure 5: The dependence of asymptotic separation between the flavour branes in units of 



R-D^/VUo with ut and magnetic field in the high temperature limit. Figure 5(a) shows the de- 



pendence of L with yx for different magnetic field strength; blue (bottom most) curve corresponds 
to H = 0, red (middle) corresponds to H = 0.5, green (top most) corresponds to H = 0.7. Figure 
|5(b)| shows the behaviour of L with applied magnetic field for different values of yx', blue (top most 
from right) curve corresponds to yx = 0, red (middle from right) corresponds to yx = 0.7, green 
(bottom most from right) corresponds to yx = 0.9. We have set d = 1. 



phase). In the presence of finite temperature there will be a first order transition to chiral 
symmetry restored phase. We will study this transition later in the next section. 




H 



0.5 1 1.5 2 



Figure 6: The dependence of constituent quark mass (measured in units of {2ira') 1 ) on the 
external magnetic field in the high temperature phase. We have fixed Ut = 0.3 and Rda = 1- 



Now the trivial solution of equation (TISI) . which is given by r' = has a different physical 
meaning from the point of chiral symmetry breaking. Since the {£, u} submanifold has the 
cigar geometry, the solutions r = ±L/2 can remain disjoint and end at u = Ut separately. 
Therefore the trivial solutions in the high temperature case preserve the full U(Nf)z x 
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U(Nf)n symmetry by remaining disjoint. In order to determine the true minimum energy 
configuration we need to compare the energies of the curved and the straight branes. We 
pursue this in the next section. 



3.2.1 The Probe Brane Profile and Chiral Symmetry Breaking 



To determine the true vacuum we consider the difference between the energies of the curved 
and straight branes, which is given by 



AS 



S 



curved 



-s t 



straight 



dyy (y 3 + HW) 1/2 
- f dyy (y 3 + HWf 2 



1 



/(l)(l+ffg d 3) 5 
f(y)(y 3 +H^d^y 



1/2 



1 



(20) 



Here AS < would mean chiral symmetry breaking, AS* > would mean chiral symmetry 
restoration and AS" = would characterize a transition from symmetry broken phase to a 
symmetry restored one. We employ numerical analysis to study this. 

It is known from, e.g., ref. [9] that for high temperature and zero magnetic field there exists 
a critical temperature beyond which the straight branes are energetically favoured, implying 
that in the dual gauge theory chiral symmetry is restored. Below the temperature, however, 
the symmetry is broken by energetically favoured curved brane pair that join together. The 



plots for the energy difference AS is shown in figures 7(a) and 7(b) for zero and non-zero 
values of magnetic field respectively. 




(a) 



(b) 



Figure 7: The dependence of zeros of AS on the magnetic field. Figure 7(a) shows the result for 



zero magnetic field and figure [7(b)] shows how the zero changes as we fix H = 1, 3, 5, 9 from above 
to below respectively. We have set d = 1. 



It is evident that the first order phase transition from chiral symmetry broken to the 
symmetry restored phase persists in presence of external magnetic field. From the zero of 
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AS we can find out the critical value of yx for which the symmetry restoration occurs. 
Now to represent the phase diagram in terms of physical quantities, we recall that there is a 
length scale L corresponding to the separation of the brane-anti-brane pair at the boundary 
So we express the chiral symmetry restoring temperature in units of 1/L using the critical 
value of ut in equation (}T9]) . The resulting phase diagram is shown in figure 

T 
0.24 

0.22 *S Restored 

0.20 




*SB 



6 8 10 12 H 



Figure 8: The phase diagram between applied magnetic field and the chiral symmetry restoring 
temperature T x sb in units of where L is the asymptotic separation between the branes. We 
have set d = 1. 



It is interesting to note that the presence of magnetic field increases the symmetry restor- 
ing temperature. In other words it promoted the spontaneous breaking of chiral symmetry. 
This fits with the general expectations from field theory (see e.g. refs.[30]) and the su- 
pergravity/probe intuition that introducing a magnetic field places more energy into the 
system; therefore in order to minimize the energy, condensates are formed (the branes bend 
more) resulting in more readily broken chiral symmetry. (It should be noted, however, that 
in this specific holographic model the identification of a quark condensate is a rather subtle 
issue (see e.g., ref. [27]).) 

We can extract some more information about the transition by studying certain thermody- 
namic quantities at the phase transition. To that end, let us note that the first order phase 
transition is accompanied by entropy density that jumps at T = T c yielding a non-zero 
latent heat as reported in ref. [IT] ; also a change in magnetization 

A_ 1 (9 (•S'curved ^straight) y~i rp \ 

S = ^ , ^latent ~ 1 c^ s j 

V R3 dT 

1 (•S'curved ^straight) /oi \ 

fi = ~v^ m • 1 ] 

The absolute free energy and any thermodynamic quantity obtained from it (such as the 
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absolute magnetization) for the two classes of embeddings (the straight and the curved 
branes respectively) are formally divergent quantities. Hence we compute the relative quan- 
tities which are finite. We studied the dependence of the change in entropy density and 



the relative magnetization numerically, and the results are shown in figures 9(a) and 9(b) 



respectively. The relative magnetization also shows a similar saturation behaviour for high 
enough magnetic field. The straight branes correspond to the melted phase where quarks 
are free whereas the curved branes correspond to the mesonic phase where quarks exists 
in the form of bound states or chiral condensates. Therefore it is expected that the chiral 
symmetry restored phase (corresponding to the straight branes) is more ionized than the 
chiral symmetry broken phase (corresponding to the curved branes). This is in accord with 



our observation that the relative magnetization is negative in figure 9(b) 
As 




H 




10 12 



H 



(a) 



(b) 



Figure 9: The behaviour of jump in entropy and change in magnetisation at the critical tempera- 

7/2 3/2 

ture where chiral symmetry is being restored. The vertical axis is evaluated in units of CU 
and we have set d = 1. 



4 Probing with Dp-Brane 

We now consider general Dp/ Dp-brane as flavours in the colour D4-brane background. The 
low temperature background is given by equation ([TJ and the high temperature is given by 
equation Also we choose to parametrize the S 4 as follows 

dVtl = del + sin 2 Q x dQ\ + sin 2 9 X sin 2 2 (d9j + sin 2 9 z d9f) . (22) 

There are several ways to place the probe brane. We list these possibilities in the table 
below (being somewhat cavalier with the use of the term QCD), as done in ref. [13] 
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To introduce flavour brane in this set up we need to introduce probe branes that extend 
in the u direction all the way up to u — oo. In general we consider introducing Dp — Dp 
brane which intersect the colour ZM-branes at r = ±L/2 at the boundary, where the 
=F sign corresponds to position of the Dp and Dp-branes respectively. It can be noted 
that the coordinate r ranges from — ttR to +irR, and therefore L G [— nR/2, +7rR/2]. By 
construction the flavour symmetry grourjfl is U{Nf)i x U (N/)r. Depending on the dimension 
of the probe brane, the number of common directions between the colour and the flavour 
branes is determined, which corresponds to the dimension of the dual gauge theory. Such 
construction has previously been discussed in ref. p2]. We follow the same approach as 
before to introduce magnetic field. By construction we can see that the magnetic field 
which has support along the (x 2 ,x 3 ) direction will have no effect on QCD 2 . 

4.1 Low Temperature Phase 

We use similar ansatz as in equation ([5]) for Dp / ' Dp-br&ne, however depending on the di- 
mension of the probe brane we now place it at 9i — 7r/2 as required following the table. 
We also denote the profile of the probe brane by r(u) (we change the choice of function 
representing the brane profile, earlier we considered to parametrize them as u(t)). The 
probe brane DBI actiorjj can be written in the following generic manner 

S Dp = ^J due-* {gu9uu9 n x - 3 (9l x + H 2 )(detS p ^ n )) 1/2 , (23) 

where g a b are the induced metric components on the world-volume of probe brane; n is the 
dimension of the gauge theory and (p+ 1) is the dimension of the probe Dp-br&ne. T is the 
temperature of the background and M is given by 

Af = /jpVRnVsp-n . (24) 

3 For the (2 + l)-dimensional case, chiral fermions do not exist. Therefore the global symmetry is just an 
U(Nf) x U(Nf) which we continue to refer as the chiral symmetry. 

4 Once again it is straightforward to verify that there will be no Wess-Zumino term in any of these cases. 
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Similar to the D8-brane case it is trivial to see that this action is independent of t(u) and 
contains only the first derivative of t(u) with respect to u. Therefore the corresponding first 
integral of motion will be the momentum corresponding to t(u), which is: 

3 r(2n-p) 



u p - 



const. 

\RdiJ f(u) J 

We see that this equation has the general structure as the L>8-brane considered in equa- 
tion only difference appears in the power of an overall factor of u, which now contains 
the general information of the dimension of the probe brane and the gauge theory. Therefore 
magnetic field would have similar effects on the profile by bending them more, promoting 
spontaneous chiral symmetry breaking. 

The action can be computed by substituting r' from equation (f25|) into equation (1231) . 
The result is 

Afdlu? +1 - n [°°, y^(l + H 2 (d/y) 3 )y-V 2 



S d p = ~ »,„ r / dy 



9sTd~^ Ji ( y ep m (1 + m (d/yf) - /(!)(! + H*d*j) 



1/2 



where we have used the rescaled variables as denoted in equation (ITU1) . The asymptotic 
separation of the brane-anti-brane system is again given by L/2 = dyr' . 

The physics of chiral symmetry breaking remains the same. As before, in the confined 
phase chiral symmetry is always broken which is enforced by the cigar geometry in {r, u] 
submanifold of the background geometry. Next we consider the high temperature case and 
study the effect of magnetic field for a general probe brane. 

4.2 High Temperature Phase 

We carry out exactly similar analysis as before. We now find two distinct class of solutions, 
namely the straight branes and the curved branes. The energy of the curved branes is given 
by 

^curved = 2NU* +1 - n d^ f 1 ^_ Z-h-^ (1 + H 2 d?z) X \ 



' ^(l + ff%)(l-^)-(l + W)(l-^)) 1/2 ' 

(25) 

The the energy of the straight branes is given by 



^straight = 2MUZ +l - n d^ % t ^^— m — 

Dp g s Td^ 2n ~^ Jo 



z 1 ' 2 . (26) 
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The true vacuum is therefore determined by minimizing AS ~ (^ u p rvod - ^ aight ); AS < 



means chiral symmetry breaking whereas AS > means chiral symmetry restoration. By 
looking at equation (|25p and ([2"oT) one can conclude that the phase diagram will depend on 
the combination of (p + 2n) for an n-dimensional gauge theory with flavour branes wrapping 
an internal (p — n) directions. Specifically speaking for {n = 4, p = 4} and {n = 3, p = 6} 
we should have exactly the same phase diagram. It can be seen that the zeros of AS depends 
on H in exactly similar way for all these examples, so the basic nature of the phase diagram 
would be universal. All there is left for the exponent (p + 2n) then is to scale the symmetry 
restoring temperature. Numerical results are shown below. 



0.30 



0.25 



0.20 



0.15 




(a) 



(b) 



Figure 10: The phase diagram for n = 4 and n = 3. Figure 10(a) shows the phase diagram for the 



4-dimensional gauge theory, the curves correspond to p = 8, 6, 4 from below to above respectively. 



Figure 10(b) shows the corresponding phase diagram for 3-dimensional gauge theory, where the 
curves correspond to p = 6, 4 from below above respectively. We have set d = 1. 



Indeed we observe a generic nature of phase diagram in such holographic models. We 
also find that the upper bound to symmetry restoring temperature that can be reached 
introducing external magnetic field depends on the dimension of the gauge theory and the 
probe brane. We can observe that for a given dimension of the gauge theory the lesser the 
number of directions wrapped by the probe brane along the internal S* 4 , the higher is the 
symmetry restoring temperature. Thus the information of the dimensions wrapped by the 
probe brane is also encoded in the phase structure. 

5 A Note on Background Electric Field 

For this section we again consider the 8-branes as flavours. We can introduce a background 
electric field by considering the following form of world-volume gauge field [29, 26J: 

A x {t,u) = (-Et + h(u)) . (27) 



19 



This means we have a non-zero constant electric field along x 1 . The function h(u) encodes 
the response of the fundamental flavours to the external field, namely it encodes the infor- 
mation of the non-zero current when flavours are free to move and therefore conduct. Note 
that here we do not introduce any chemical potential, therefore in the gauge theory there is 
no a priori candidate for carrying the charge. However, there could still be current caused 
by pair creation in presence of the electric field. We comment on some observations and 
expectations henceforth. 

For simplicity we assume 27ra f = 1. We choose the same ansatz for the probe brane profile 
as equation ([5]). The DBI action can be computed to be 

S D8 = cfdudte-«G„(tetstf[G tt G xx9u , + G u h°- gu ,E>r, 
9uu — G uu + t G TT , (28) 

where C = (Nf^/CV^Vs^), G^ v is the background metric and g uu is the induced metric 
component on the world-volume of the D8-brane along w-direction. 

Now we will have two constants of motion corresponding to the function r(u) and h(u) as 
follows 

Ce-*y/MS A G xx {GttG xx - E 2 )G TT r' fi 
yjG u h n - (E 2 - G tt G xx )(G uu + G TT r' 2 ) 



Ce-*y/Ms:G xx G t th' = j 

y/Guh* - {E 2 - G tt G xx ){G uu + G TT r' 2 ) 

The constants B and J are related to the minimum radial distance along the world-volume 
of the probe brane (Uq as before) and the gauge theory current (a similar identification 
made in ref. [29J holds here also). 

Solving equation (1291) numerically to look for all possible solutions is a difficult problem. 
Nevertheless, it can be shown from equation fl29|) that for confined phase there can be 
only the joined solution (curved as before), for which the current vanishes if we impose the 
condition that the brane-anti-brane pair join smoothly at Uq. This can be shown by formally 
solving equation fl29|) to obtain r{u) and h(u) as a function of u and the constants B and 
J. If we expand the solution around the joining point u = Uq and demand that r' — > oo 
as u — > Uq, we obtain that the current identically vanishes. However, if the branes do join 
but not smoothly, this is no longer necessarily true. One could, in this case, have cusp-like 
solutions for which the brane-anti-brane pair join at an angle and in order to stabilize the 
system it is necessary to consider the inclusion of a bunch of fundamental strings extending 
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from the joining point u = Uq to u = Ukk- This leads us to a construction much like 
in refs. [T4"l [T5| [16] , where the effect of baryons was considered. Therefore the symmetry 
broken phase may have a non-zero current carried by the baryons. The chiral symmetry is 
always broken in this phase forced by the topology of the background. 

For the deconfined phase, it can also be shown in a similar way that there exists curved 
solutions joining smoothly at some Uq, which has zero current modulo the caveat mentioned 
in the last paragraph. This is intuitive from the gauge theory point of view, since we are 
in a chiral symmetry broken phase therefore there is no charge carrier present to conduct 
(ignoring the possibility of a baryon current). The possible effect of pair creation is dimin- 
ished by the existence of quark bound states in the symmetry broken phase. However, for 
straight branes we expect non-zero current to flow. 

A familiar fact from studying flavours in electric field tells us that the presence of electric 
field induces a so called "vanishing locus" for the probe DBI action. The "healing" procedure 
(e.g., in ref. [26]) is to give a non-zero vev to the current. In practice this is obtained by 
substituting the functions t(u) and h(u) from equation (|29|) in favour of the constants B 
and J in the action in equation (1281) and demanding the reality condition for the action for 
u G [Ut, oo]. The condition leads to the following two equations 

(G U G XX — E 2 ) = , 

B 2 e 2 *G tt - (G tt G xx - E 2 ) G TT {C 2 {detS A )G tt Gl x - e^J 2 ) = . (30) 

The two expressions are the terms in the on-shell action that go to zero in the numerator 
and the denominator respectively. It can be shown from equation f[3"Uj) that in order to have 
J 7^ one has to have B = 0, which corresponds to the straight brane solutions. 

For a given electric field, we can determine the position of the vanishing locus u e h from 
the first condition in equation (I3"0"j) . Knowing u e % we can then extract the current J, and 
therefore the conductivity using J = oE. To express the conductivity in terms of physical 
quantities let us recall that 

fi 8 = (2tt)-V-2 , T 



Now we will restore the factors of (27ra') and also set R? D4: = nXa' = 1. Combining everything 
we get the following expression 

4 , / 27 E \ 1/3 

o^AJW 2 (1 + 3^5) ■ (3D 

This conductivity is due to the melting of mesons at high temperature and pair creation 
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mediated by the electric field. Since we get this formula using DBI action, it captures 
non-linear behaviour of the conductivity with respect to the electric field. 

The effect of the electric field on chiral symmetry breaking should be to reduce the sym- 
metry restoring temperature by polarizing the bound states into constituent quarks. The 
analysis above points to the fact that electric field works as expected from field theory per- 
spective. However, the energy consideration does not lead to the expected result, because 
as electric field increases u e h also increases, therefore the straight branes which extend all 
the way down to Ut has more DBI action energy as compared to their curved counterparts, 
which can only extend down to u e h- We hope to address this issue in future. 

6 Conclusion 

We have extended the study of Sakai-Sugimoto model to include the presence of external 
electric and magnetic fields, examining the dynamics of the flavour sector in (an analogue 
of) the "quenched" approximation. We have seen that external magnetic field helps in chiral 
symmetry breaking. This particular effect of an external magnetic field has been referred 
to as magnetic catalysis in field theory literature (see e.g., refs. [30J ) . Our results and 
observations are consistent with results from those approaches. We found that the chiral 
symmetry restoring temperature increases with increasing magnetic field. We have further 
observed that such holographic models have an upper bound for the symmetry restoring 
temperature depending on the dimension of the gauge theory and the probe brane. 

We briefly studied the effect of external electric field, and though we have not explored all 
of the details, we expect the dynamics to be also consistent with the field theory intuition, 
although (as we did for magnetic field here) the precise details should be interesting to 
uncover. The presence of non-zero electrical conductivity is also an avenue of further study. 
It would also be interesting to study Hall effect in such models when both electric and 
magnetic fields are present. 

Another important avenue would be to study the meson spectra in the presence of these 
external fields, where we expect to see effects such as Zeeman splitting, the Stark effect and 
so forth (as observed in the D3/D7 model in some of refs. [2S])- We leave these and many 
other interesting aspects for future exploration. 



22 



Acknowledgments 



We would like to thank Tameem Albash, Nikolay Bobev, Veselin Filev and Ramakrishnan 
Iyer for discussions. This research was supported by the US Department of Energy. 



References 

[1] Juan M. Maldacena, "The large N limit of superconformal field theories and 



supergravity,Mdu. Theor. Math. Phys. 2 (1998) 231-252, hep-th/9711200 



[2] E. Witten, "Anti-de Sitter space and holography," Adv. Theor. Math. Phys. 2 (1998) 
253-291, |hep-th/9802150| 

[3] S. S. Gubser, I. R. Klebanov and A. M. Polyakov, Phys. Lett. B 428, 105 (1998) 
[arXiv:hep-th/9802109| . 

[4] O. Aharony, S.S. Gubser, J. Maldacena, H. Ooguri and Y. Oz "Large N field theories, 



string theory and gravity," Phys. Rept. 323 (2000) 183-386, hep-th/9905111 



[5] E. Witten, "Anti-de Sitter space, thermal phase transition, and confinement in gauge 
theories," Adv. Theor. Math. Phys. 2 (1998) 505-532, |hep-th/980313lj 

[6] A. Karch and E. Katz, "Adding flavor to AdS/CFT," JEEP 06 (2002) 043, 



hep-th/ 0205236. 



[7] T. Sakai and S. Sugimoto, "Low energy hadron physics in holographic QCD," Prog. 



Theor. Phys. 113 (2005) 843-882, hep-th/0412141 



T. Sakai and S. Sugimoto, "More on a holographic dual of QCD," Prog. Theor. 



Phys. 114 (2006) 1083-1118, |hep-th/0507073 



[8] E. Antonyan, J. A. Harvey, S. Jensen and D. Kutasov, "NJL and QCD from string 
theory," |hep-th/06040T7[ 

E. Antonyan, J. A. Harvey and D. Kutasov, "The Gross- Neveu model from string 
theory," |hep-th/0608T49| 



E. Antonyan, J. A. Harvey and D. Kutasov, "Chiral symmetry breaking from 
intersecting D-branes," hep-th/0608177| 

[9] O. Aharony, J. Sonnenschein and S. Yankielowicz, "A holographic model of 
deconfinement and chiral symmetry restoration," |hep-th / 0604161, 



23 



[10] R. Casero, A. Paredes and J. Sonnenschein, "Fundamental matter, meson 
spectroscopy and non-critical string / gauge duality," JHEP 0601, 127 (2006) 
|arXiv:hep-th/0510lT0) . 

K. Peeters, J. Sonnenschein and M. Zamaklar, "Holographic melting and related 
properties of mesons in a quark gluon plasma," Phys. Rev. D 74, 106008 (2006) 
|arXiv:hep-th/0606195] . 

K. Peeters, J. Sonnenschein and M. Zamaklar, "Holographic decays of large-spin 
mesons," JHEP 0602, 009 (2006) |arXiv:hep-th/0511044| . 

[11] A. Parnachev and D. A. Sahakyan, "Chiral phase transition from string theory," 
Phys. Rev. Lett. 97, 111601 (2006), |hep-th/0604173l 

A. Parnachev and D. A. Sahakyan, "Photoemission with chemical potential from 
QCD gravity dual," Nucl. Phys. B 768, 177 (2007) |arX iv:hep-th/0610247l . 

[12] N. Horigome and Y. Tanii, "Holographic chiral phase transition with chemical 
potential," JHEP 01, 072 (2007), jhep-th/0608198 

[13] D. Gepner and Shesansu Sekahr Pal, "Chiral symmetry breaking and restoration from 



holography," hep-th/0608229 



[14] M. Rozali, H. H. Shieh, M. Van Raamsdonk and J. Wu, "Cold Nuclear Matter In 
Holographic QCD," JHEP 0801, 053 (2008) [arXiv: 0708 .13221 [hep-th]]. 

[15] O. Bergman, G. Lifschytz and M. Lippert, "Holographic Nuclear Physics," JHEP 
0711, 056 (2007) [ arXiv:0708.0326l [hep-th]]. 

[16] J. L. Davis, M. Gutperle, P. Kraus and I. Sachs, "Stringy NJL and Gross-Neveu 

models at finite density and temperature," JHEP 0710, 049 (2007) [arXiv:0708.0"589l 
[hep-th]]. 

[17] O. Aharony, K. Peeters, J. Sonnenschein and M. Zamaklar, "Rho meson condensation 
at finite isospin chemical potential in a holographic model for QCD," JHEP 0802, 
071 (2008) [ arXiv:0709.3948l [hep-th]]. 

[18] A. Parnachev, "Holographic QCD with Isospin Chemical Potential," JHEP 0802, 062 
(2008) [ arXiv:0708.3T70l [hep-th]]. 

[19] M. Kruczenski, D. Mateos, R. C. Myers, and D. J. Winters, "Towards a holographic 



dual of large-N(c) QCD," JHEP 05 (2004) 041, hep-th/0311270 



24 



[20] J. Babington, J. Erdmenger, N. J. Evans, Z. Guralnik, and I. Kirsch, "Chiral 

symmetry breaking and pions in non-supersymmetric gauge / gravity duals," Phys. 
Rev. D69 (2004) 066007, |hep-th/03060l8l 

[21] M. Kruczenski, D. Mateos, R. C. Myers, and D. J. Winters, "Meson spectroscopy in 



AdS/CFT with flavour," JEEP 07 (2003) 049, hep-th/0304032 



[22] T. Albash, V. Filev, CV. Johnson and A. Kundu, "A topology-changing phase 
transition and the dynamics of flavour," hep-th/0605088, 

[23] D. Mateos, R. C. Myers, and R. M. Thomson, "Holographic Phase Transitions with 
Fundamental Matter," |hep-th/0605046 



[24] A. Karch and A. O'Bannon, "Chiral transition of N = 4 super Yang-Mills with flavor 
on a 3-sphere," |hep-th/0605T20l 

[25] K. Ghoroku, T. Sakaguchi, N. Uekusa, and M. Yahiro, "Flavor quark at high 
temperature from a holographic model," Phys. Rev. D71 (2005) 106002, 



hep-th/0502088 



[26] V. G. Filev, C. V. Johnson, R. C. Rashkov and K. S. Viswanathan, "Flavoured 
large N gauge theory in an external mag netic field," |hep-th/ 0701001, 
T. Albash, V. G. Filev, C. V. Johnson and A. Kundu, "Finite Temperature Large N 
Gauge Theory with Quarks in an External Magnetic Field," larXiv:0709. 15471 [hep-th]. 
T. Albash, V. G. Filev, C. V. Johnson and A. Kundu, "Quarks in an External Electric 
Field in Finite Temperature Large N Gauge Theory," larXiv:0709. 15541 [hep-th]. 
J. Erdmenger, R. Meyer and J. P. Shock, "AdS/CFT with Flavour in Electric and 
Magnetic Kalb-Ramond Fields," JHEP 0712, 091 (2007) jarXiv:0709.155TI [hep-th]]. 
S. Penati, M. Pirrone and C. Ratti, "Mesons in marginally deformed AdS/CFT," 
JHEP 0804, 037 (2008) |arXiv:0710.4292l [hep-th]]. 

[27] O. Bergman, S. Seki and J. Sonnenschein, "Quark mass and condensate in HQCD," 
JHEP 0712, 037 (2007) [ arXiv: 0708 .28391 [hep-th]]. 

[28] O. Bergman, G. Lifschytz and M. Lippert, "Response of Holographic QCD to Electric 
and Magnetic Fields," larXiv: 0802 .37201 [hep-th]. 

A. Karch and A. O'Bannon, "Metallic AdS/CFT," JHEP 0709, 024 (2007) 
|arXiv:0705.3870l [hep-th]]. 

A. O'Bannon, "Hall Conductivity of Flavor Fields from AdS/CFT," Phys. Rev. D 76, 
086007 (2007) [ arXiv: 0708 .19941 [hep-th]]. 



25 



[30] V. A. Miransky, "Dynamics of QCD in a strong magnetic field," 
arXiv:hep-ph/0208180"l 

G. W. Semenoff, I. A. Shovkovy and L. C. R. Wijewardhana, "Universality and the 
magnetic catalysis of chiral symmetry breaking," Phys. Rev. D 60, 105024 (1999) 
|arXiv:hep-th/9905iT6| . 

V. P. Gusynin, V. A. Miransky and I. A. Shovkovy, "Dimensional reduction and 
catalysis of dynamical symmetry breaking by a magnetic field," Nucl. Phys. B 462, 
249 (1996) |arXiv:hep-ph/9509320] . 



26 



